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Abstract 



The shape of the V — I characteristics leading to a peak in the differential 
resistance = dV/dl in the vicinity of the order-disorder transition in NbSe2 
is investigated, is large when measured by dc current. However, for a 
small I ac on a dc bias r^ decreases rapidly with frequency, even at a few 
Hz, and displays a large out-of-phase signal. In contrast, the ac response 
increases with frequency in the absence of dc bias. These surprisingly opposite 
phenomena and the peak in are shown to result from a dynamic coexistence 
of two vortex matter phases rather than from the commonly assumed plastic 
depinning. 

PACS numbers: 74.60.Ec, 74.60.Ge, 74.60.Jg 
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Measurement of V — I characteristics or of the differential resistance r^ = dV/dl in 
superconductors is a very common method to investigate vortex dynamics and to identify 
various possible pinned and moving phases of the vortex matter. In particular, the specific 
shape of V — I that leads to a peak in dV/dl has attracted much attention ||T|-|TT|1. In 
numerical simulations this peak in is usually ascribed to plastic vortex depinning followed 
by dynamic ordering of the lattice |8Hill . However, magnetic decorations and SANS studies, 
which have observed the plastic deformation of the lattice near the depinning followed by 
the dynamic ordering fl2Hl4l1, did not find the predicted peak in dV/dl fl2|Jl3l| . In addition, 
in clean systems like NbSe2 the peak in dV/ dl is present surprisingly only in the lower part 
of the peak effect (PE), whereas in the rest of the H — T phase diagram the V — I curves 
are concave upward with no peak in dV/dl P,|3|,[T5|1 . Moreover, in the same region where 



the peak in r^ is observed, a number of anomalous vortex matter properties were recently 
found, including the striking observation that for an ac current the apparent vortex mobility 
increases rapidly with frequency [p^-|20|. Several ideas and models have been proposed 
in which the ac agitation facilitates the plastic vortex depinning. Yet another important 
paradox has received little attention: when the ac current is superposed on a dc bias the 
opposite behavior is observed, i.e., the apparent vortex mobility decreases with frequency, 
even at frequencies as low as several Hz || . None of the models that describe the peak in 
or the mobility enhancement with frequency fl8|- |lTl . |T6| has resolved this apparent paradox. 

In this paper we demonstrate that the shape of the V — I curves, the peak in dV/dl, 
and the opposite frequency dependencies stem, instead, from a dynamic coexistence of two 
vortex matter phases. A highly-pinned metastable disordered phase (DP), generated at the 
sample edges, anneals into a weakly-pinned equilibrium ordered phase (OP) in the bulk of 
the sample fTT| , |2"T| , p2f . The specific shape of the dc V — I curves results from the fact that 
most of the sample is in the metastable DP at low currents, but in the OP at high currents. 
This dynamic transformation results in a peak in dV/dl. If measured by a small ac current 
superposed on a dc bias, we find that the peak in dV/dl decreases with frequency and 
displays a unique out-of-phase signal due to the slow transformation process. In contrast, 



for an ac current with no dc bias only the edges of the sample are contaminated by the DP, 
resulting in the opposite behavior, i.e., the voltage response grows with the ac frequency. 

Transport measurements were carried out on several Fe-doped (200 ppm) NbSe2 crystals 
in strip-like four-probe configuration in applied field H\\ c-axis. The data presented here 
are for a crystal 2 x 0.4 x 0.04 mm 3 with T c = 5.6 K, H c2 (A.2K) = 1 T, and the PE field 
H P (4.2K) = 0.53 T. Very low contact resistance of ~ 10 mQ was achieved with large current 
contacts of Au evaporated onto a freshly cleaved surface. By immersing the crystals in 
liquid He, currents up to 100 mA could be applied with negligible heating. Square wave or 
sinusoidal I ac was used and the corresponding V ac was measured by a lock-in amplifier. 

Figure 1 shows V ac vs. I ac and V dc vs. I dc in the vicinity of the order-disorder transition 
in the lower part of the PE at H = 0.44 T. The dc curve (solid circles) starts to increase 
in concave form above 18 mA and then rapidly turns convex. At higher currents linear 
flux flow behavior is obtained. Although the form of the dc V — I curve may seem rather 
conventional, we find that in NbSe2 the convex shape is present only in the lower part of the 
PE, while in the rest of the phase diagram the curves are always concave upward, consistent 
with previous reports PJ3|,|I5||. V ac vs. I ac measured at various frequencies in Fig. 1 are 
remarkably different from the dc V — I. Even at a frequency as low as 1 Hz the apparent 
J c is much lower and the voltage response below 20 mA is strongly enhanced. Furthermore, 
the apparent vortex mobility increases rapidly with frequency, as noted previously [p!^ -[T8|. 

Figure 2 shows the differential resistance r d measured by superimposing a small I ac (0.1 
to 1 mA) on I dc , r d = V ac / I ac , along with r d = dV/dl obtained by numerical differentiation 
of the dc V - I at H = 0.31 T. At this slightly lower field within the PE the dc V - 
I characteristic (solid curve) is more gradual and turns convex above ~ 20 mA. At the 
inflection point, dV/dl displays a large peak reaching three times the flux flow ry. The 
remarkable result here is that the ac r d is significantly different from the dc value || and it 
decreases rapidly with ac frequency. Note that the peak in ac r d is suppressed to about half 
of the dc value already at / as low as 3 Hz, indicating the existence of very long characteristic 
timescales that are even longer than the vortex transit time across the sample, as described 
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below. None of the microscopic bulk mechanisms, such as plastic depinning or dynamic 
ordering, can account for such long timescales. Moreover, due to the high sensitivity of lock- 
in technique, the ac vs. Id c is commonly integrated to derive the full V — I characteristics 
@||. The important conclusion from Fig. 2 is that the integral of the ac is not equal to 
the dc V - I. 

Figures 1 and 2 display a striking qualitative difference: for a pure I ac the measured V ac 
increases with frequency, whereas for a small I ac superposed on a larger the resulting 
V ac decreases with /. The lower panel of Fig. 2 shows another puzzling aspect of the data, 
which is a large out-of-phase component that appears only in the non-linear part of the 
V — I and is absent in the flux-flow region. An out-of-phase signal is a common feature in 
ac susceptibility measurements, where the amplitude of the current induced in the sample 
and the dissipation level depend on the excitation frequency. In transport measurements, 
in contrast, the amplitude of the current is fixed by the external circuitry and therefore the 
voltage response, as a rule, is frequency independent and does not show any out-of-phase 
signal at low frequencies. To the best of our knowledge this is the first published report of 
an imaginary r<£, which further emphasizes the anomalous vortex dynamics in the lower part 
of the PE. Note that the out-of-phase is non-monotonic with frequency: it vanishes in 
the limit of high and low / and is largest for the 22 Hz data. 

We now discuss the results in view of the recent understanding that the PE reflects a 
disorder- driven first-order phase transition from a weakly pinned OP (Bragg glass) with 
a low critical current density J° rd into a strongly pinned DP with a high J™* 8 p^-p5|. 
Local measurements have demonstrated that below the transition, in the lower part of 
the PE, in the presence of transport current a supercooled metastable DP is formed at 
the sample edge because of non-uniform vortex penetration through the surface barriers 
p!7| , |2T|j22]] . As the vortex lattice moves across the sample, the metastable DP with a high 
concentration of dislocations gradually anneals into the dislocation-free OP. We describe 
the annealing stage of the DP by its local critical current density J c (x) which has a non- 
equilibrium excess value J c {x) = J c {x) — J° rd relative to the fully annealed OP. Since in low- 
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temperature superconductors thermal activation is negligible, the sole annealing mechanism 
of the metastable DP is through a current-driven displacement that allows rearrangement 
and disentanglement of the vortices during the motion. We therefore assume, for simplicity, 
that the relative annealing of J c upon displacement by a small Ax is given by Ax/L r , where 
L r is a characteristic relaxation length over which the DP anneals into the OP. Since the 
lattice flows with velocity v, J c at x + Ax and at time t + At = t + Ax/v is thus described by 
J c {x + Ax, t + Ax/v) = J c (x, t)(l — Ax/L r ), which leads to the partial differential equation 
of the annealing process 

dJ c (x, t)/dx + {l/v)dJ c (x, t)/dt = -J c {x, t)/L r (v), (1) 

with a boundary condition at x = 0, where vortices penetrate into the sample, of J c (0, t) = 
jdis _ jord^ aS p ec t Q f ^he annealing process is that the relaxation length L r crucially 

depends on the displacement velocity v. Fast transient measurements [26[] show that at 



low velocities L r is large, whereas for a potential landscape that is strongly tilted by a 
large driving force the disentanglement is very rapid, so that empirically L r ~ L (vq/v) v = 
Lo(yo/V) v . Here r\ is typically in the range of 1 to 3, L , v , and V are scaling parameters, 
V = vBl is the measured voltage drop, B is the magnetic field, and / is the distance between 
the voltage contacts. We now demonstrate that these simple assumptions describe all the 
essential experimental observations. 

We first analyze the time independent behavior. The dc solution of Eq. |l| is Jc°{x) = 
(^jdis _ j° rd ^ exp(—x/L r ) + J° rd , as shown schematically in Fig. 3 inset. By integrating over 
the width W we obtain the total critical current of the sample I c = d J dc (x)dx: 

I c {L r ) = (jf s -J: rd ){l-e~^(vl}L r (V)d + r c rd , (2) 

where I° rd = J° rd Wd and d is the sample thickness. Note that I c depends on L r , which in 
turn depends on voltage V. This property is central to the described phenomena: Coexis- 
tence of the DP and OP results in an inhomogeneous sample and moreover the degree of the 
inhomogeneity, J c (x), changes with vortex velocity. As a result, the total I c of the sample 
is not fixed, but rather changes with voltage. 



The dc V — I characteristics can be derived as following. We write for simplicity the 

V — I of the OP as V = rj(I — I° rd ), where rf is the flux-flow resistance. Similarly, when the 
entire sample is in the DP, V = r/(J — I dls ) with l dts = J diS Wd. These two asymptotic V — I 
solutions are shown by the dashed lines in Fig. 3. Here we have assumed for simplicity that 
the flux-flow resistance ry is the same for any of the phases. As a result, when the two phases 
coexist V = rf(I — I c ), where I c is given by Eq. § and is voltage dependent through L r (V). 
Hence an analytical I(V) relation can be written directly as I = V/rj + I c (L r (V)), and the 
resulting non-linear V — I characteristic is shown by the solid curve in Fig. 3. At very low 
voltages L r is larger than the sample width, namely the entire sample is contaminated by the 
DP, and hence the V — I initially follows the asymptotic dashed line of the DP with J c = I dls . 
At high vortex velocities L r becomes very short, most of the sample is in the OP and the 

V — I approaches the asymptotic line of the OP with I c = I° rd . In the crossover region 
a specific shape of the curve with an inflection point is obtained alike the experimental dc 

V — I curves in Figs. 1 and 2. This shape is the result of a continuous decrease of the total 
I c of the sample from J c = I dls to J c = I° rd with increasing current. The exact curvature 
in the crossover region depends on the parameters T], L , and Vq (see Fig. 3 caption) and 
may be either gradual or very steep, and may even obtain a negative slope resulting in an 
S-shaped characteristic as found recently in the lower part of the PE \ ^% . 

In order to understand the frequency dependence of the differential resistance shown 
in Fig. 2 and in more detail in Fig. 4a, we solve Eq. [I] for a small periodic veloc- 
ity perturbation v = Vd c + v ac e luJt caused by applied current l dc + I ac e luJt - In this case 
L r {y) = L dc + (dL r /dv)v ac e luJt . Following a simple calculation and keeping only the 
linear terms in v ac , we obtain J c (x,t) = J dc (x) + J^ c (x)e tuJt , where Jc C {%) = {Jc ls ~ 
J° rd )(dL r /dv)(v dc v ac / L 2 dc )e~ x/Ldc (l - e~ iujx/vdc )/iuj. Note that J c ac (x) is negative because 
of a negative dL r /dv, reflecting the fact that the system becomes more ordered with increas- 
ing v. In the limit of low frequency, J c (x,t) slowly varies between two extreme dc solutions 
determined by L dc ± (dL r /dv)v ac , as shown in Fig. 3 inset. At the minimum value of the 
current, — I ac , the vortex velocity is lowest, L r is largest, and hence J c {x) is highest. At 
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the maximum of the current, I dc + I ac , L r is smallest and the sample is 'cleanest', resulting in 
a large enhancement in vortex velocity. Hence the large r d at low / arises from the varying 
contamination of the sample which significantly amplifies the voltage response. However, 
modification of the sample contamination is a slow process because the only mechanism of 
enhancement of the local disorder in the bulk is by transporting a more disordered lattice 
from the edge of the sample. This results in characteristic timescales comparable with the 
vortex transit time across the sample r t = W/v (or even longer, see below). 

To obtain the full frequency dependence of r d we calculate = d J^ c (x)dx, and 
by noting that V ac = r f (I ac - J c ac ), we find r d = V ac /I ac = (l/r f + I* c /V ac y l . This result 
shows that the enhancement of r d results from the fact that J" c is negative, and depending 
on the L r (v) parameters, r d can become infinite and even negative for the case of S-shaped 
V — I curves. The full analytical expression for r d is too extensive to be presented here. 
Figure 4b therefore shows the calculated real and imaginary parts of r d /rf vs. frequency at 
the operating point in Fig. 3. At low frequencies the in-phase r d is large and it decreases 
towards rj when / approaches the transit frequency f t = l/r t . We can understand this 
behavior as follows. As indicated by the arrows in Fig. 3, at low frequencies r d is given 
by dV/dl which can be very large depending on the specific shape of the dc V — I curve. 
However, the surprising result here is that at higher frequencies, in contrast to the common 
belief, the experimental r d = V ac / I ac does not measure the true dV/ dl. At high / the local 
degree of disorder J c (x, t) cannot adjust to the rapid variations and therefore J c (x, t) is fixed 
at J^ c (x), and J^ c {x) and I® c vanish. As a result, at high frequency the ac signal, instead 
of following the dc curve, follows a trajectory shown by the dotted line in Fig. 3. This line 
is the V — I characteristic of a sample with a fixed I c = If 1 , V = r/(7 — If), resulting in 
r d = rf, as indeed observed experimentally in Figs. 2 and 4a. 

The rearrangement of the bulk disorder always lags the external ac drive, giving rise 
to a pronounced imaginary r d component (Fig. 4b) that has a maximum at intermediate 
frequencies. The qualitative agreement between the results of our simplified calculations in 
Fig. 4b and the experimental data in Fig. 4a is remarkable, including the frequency scale: 
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The experimental transit frequency ft = Vdc/BWl, marked by the arrow in Fig. 4a, was 
obtained by a direct measurement of Vd c , and the maximum of the out-of-phase r d occurs 
at about 0.1ft, in excellent agreement with the calculated behavior 



Finally, we comment briefly on the response to a large I ac with no Idc as shown in Fig. 
1. In contrast to the previous situation, here all the metastable DP exits and reenters 
the sample during every ac cycle, as shown previously ||17|| . The detailed analysis of this 
case is much more complicated since J c {x, t) changes significantly during the ac cycle and 
is different for the entering and exiting edges, resulting in highly nonlinear behavior that 
cannot be treated analytically. Qualitatively, however, the maximum depth to which the DP 
penetrates the sample from each edge during the corresponding half period of the ac cycle 
is x a d c = v/2f, where v is the average vortex velocity during the half cycle, while the central 
part of the sample remains ordered. Therefore the fraction of the sample occupied by the 
DP decreases with / as x a d c /W oc 1/ f. Consequently, the volume of the OP increases with 
frequency, the integrated I c decreases, and V ac increases. At sufficiently high /, practically 
the entire sample becomes ordered so that the experimental V ac — I ac in Fig. 1 follows the 
asymptotic dashed OP curve in Fig. 3. Thus a high-frequency measurement of V ac — I ac is a 
useful method to reduce edge contamination and to approximate the true V — I of the OP. 

In summary, the process of edge contamination by the metastable disordered phase in 
the vicinity of the order- disorder transition is shown to explain the convex shape of the 
V — I characteristics in the lower part of the PE, the large difference between V^ c — Id c arid 
V ac — I ac curves that grows with frequency, the peak in the differential resistance and its 
rapid suppression with frequency, the out-of-phase signal in ac transport measurements, and 
the opposite frequency dependence in the different cases of ac drive. 
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FIGURE CAPTIONS 

Fig. I. V ac — lac characteristics at various frequencies and V dc — I dc (•). The voltage 
response increases with ac frequency. 

Fig. 2. In-phase and out-of-phase differential resistance r d at dc (•) and various ac 
frequencies (left axis), and dc V — I characteristic (right axis). In-phase r d decreases with 
frequency, while the out-of-phase r d is maximal at intermediate frequencies. 

Fig. 3. Theoretical dc V - I characteristic (solid line) with I°J d = 2 mA, If 8 = 18 mA, 
r f = 2 mQ, W = 400 fim, and L r (V) = L^Vo/Vf with v = 2, L = 200 fim, and V = 30 
/iV. The frequency dependence of the differential resistance at the operating point (•) is 
shown in Fig. 4b. Inset: Schematic sample geometry and J c (x) across the sample width. 
J c {x) decreases with increasing current. 

Fig. 4. (a) Frequency dependence of the in-phase (•) and out-of-phase (o) r d at I dc = 22 
mA near the peak of r d in Fig. 2. (b) Calculated r d /rj vs. f/f t at the operating point in 
Fig. 3. 
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Fig. 1, Paltiel etal. 




Fig. 2, Paltiel etal. 




Fig. 4, Paltiel etal. 



